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Abstract 

A statistical field theory is developed to explore the density of states and spatial 
profile of 'tail states ' at the edge of the spectral support of a general class of disordered 
non-Hermitian operators. These states, which are identified with symmetry broken, 
instanton field configurations of the theory, are closely related to localized sub-gap 
states recently identified in disordered superconductors. By focusing separately on the 
problems of a quantum particle propagating in a random imaginary scalar potential, 
and a random imaginary vector potential, we discuss the methodology of our approach 
and the universality of the results. Finally, we address potential physical applications 
of our findings. 

1 Introduction 

Over recent years the statistical properties of stochastic non-Hermitian operators have 
come under intense scrutiny § |, |[ g | |[ |7[ 1 1 [lO[ [l^ [l], |l|, |l§ 



Operators of this kind appear in a number of physical applications ranging from random 
classical dynamics, and statistical physics, to phase breaking and relaxation in quantum 
dynamics. For example, a straightforward mapping shows the statistical mechanics of a 
repulsive polymer chain can be described in terms of the quantum dynamics of a par- 
ticle subject to a random imaginary scalar potential j2^, |2^. Similarly, the statistical 
mechanics of flux lines in a type-II superconductor pinned by a background of impuri- 
ties can be described as the quantum evolution of a particle in a disordered environment 
subject to an imaginary vector potential Finally, the diffusion of a classical particle 
advected by a random velocity field is described by a linear non-Hermitian operator, the 
"Passive Scalar" equation . More generally, the dynamics of various classical systems 
can be expressed in terms of random Fokker-Planck operators (for a review see, e.g., 
Refs. H, H). 

Non-Hermitian operators exhibit several qualitatively new phenomena which discrim- 
inate their behaviour from those of their random Hermitian counterparts. Firstly, while 
the eigenvalues of Hermitian operators are real, the spectrum of non-Hermitian operators 
is generically complex, bound by some region of support in the complex plane. The second 
striking distinction concerns the localization properties of the eigenfunctions. It is well 
established that a weak random impurity potential brings about the localization of all 
eigenstates of a Hermitian operator in dimensions of two and below By contrast, 

a constant imaginary vector potential is sufficient to delocalize states of the disordered 
system even in one-dimension! This phenomenon, which was reported and explained by 
Hatano and Nelson finds extension to higher dimensio ns p^. 

Beginning with early work on random matrix ensembles 1 26 , 2^ , a variety of techniques 



have been developed to study the statistical properties of stochastic non-Hermitian oper- 
ators 0, III, |,|, I, I, I |l^, In], |12|, [18[ Some of the techniques 
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have been based on the random matrix theory |l|, ^, ||, |[ while others 

have involved the development of perturbative schemes, such as the self-consistent Born 
approximation in the diagrammatic analysis (e.g., Ref. |^). In this paper we will be 
concerned with adapting a third approach which involves the refinement of field theoretic 
techniques which have proved to be so useful in the description of random Hermitian 
operators. 

Recently, attention has been directed towards the study of 'tail states' which persist at 
the edge of the spectral support of the non-Hermitian system. Indeed, there are already 
indications |2g| , ^ that an important role can be played by those parts of the spec- 
trum which are populated by rare states. Now, in stochastic Hermitian systems a small 
tail of states accumulate below the band edge tightly bound to rare or 'optimal fluctu- 



ations' of the random potential |31, |3^, ^ . These so-called 'Lifshitz tail' states are 



typically smooth, node-less, localized functions which inhabit regions where the potential 
is particularly shallow and smooth. By contrast, the character of the complex spectrum 
in non-Hermitian operators allows for the existence of tail states that occupy the entire 
region which bounds the support of the spectrum. 

The difference is not incidental: tail states in the non-Hermitian system can exhibit 
features characteristic of the quasi-classical states within the bulk spectrum. To emphasize 
the point, consider the spectrum of, say, a two-dimensional free quantum particle subject 
to a complex scalar potential which involves both real and imaginary components, 

H = -^+iV{r) + W{r) , (1) 
2m 

where p = —ihV. If the imaginary potential V is absent, a region of localized tail states 
accumulate below the band edge (e < 0) due to optimal configurations of the potential. 
At energies e greatly in excess of the inverse scattering time /i/r, bulk states are only 
'weakly localized' with a localization length ^loc greatly in excess of scattering mean-free 
path £ and wavelength A. This is the quasi-classical regime where mechanisms of quantum 
interference strongly affect the spectral and transport properties of the system p5[ . 

On this background, let us now suppose that a small imaginary component of the 
potential is restored. Being non-Hermitian, the eigenvalues migrate into the complex 
plane occupying a finite region of support centered on the real axis (see Fig. |l|). At the 
level of the mean-field analysis (detailed below), complex eigenvalues are confined to a 
sharp region of support in the complex plane. However, as with Lifshitz states, optimal 
fluctuations of the random potential(s) generate states which lie outside the region of 
support (see Fig. ^. What is the nature of these states? Are they locahzed and, if 
so, over what length scale? Here, in contrast to the Hermitian system, tail states must 
accommodate fast fluctuations of the wave function at the scale of the wavelength of the 
Hermitian system. 

To our knowledge, the problem of tail states in the non-Hermitian system was first 
considered in a recent work by Izyumov and Simons . Refining an approach developed 
in Refs. Q to explore tail states in Hermitian systems, the authors introduced a non- 
perturbative scheme to investigate states at the edge of the support in non-Hermitian 
systems. This instanton technique revealed the existence of tail states at the edge which 
are localized on a length scale greatly in excess of the wavelength — tail states in the 
non-Hermitian system are quasi-classical in nature. Building on the work of Efetov | |35[ |, 
the aim of the present paper is to develop a supersymmetric field theory to explore the 
spatial profile of the states in the tail region, and to determine the complex density 
of states (DoS) in the vicinity of the band edge. In doing so, we will reveal a close 
correspondence between tail states in the non-Hermitian problem and sub-gap states in 
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Figure 1: Complex eigenvalue spectrum for several realizations of the random non-Hermitian 
Hamiltonian (jl|), where the potentials V and W are both drawn from a Gaussian 5-correlated ran- 
dom impurity distribution. Here we have included a lattice of 29 x 29 q-points in the momentum 
space. 

weakly disordered superconductors. Moreover, constraints within the present scheme 
reveals the universal character of tail states in the non-Hermitian system showing that, 
after a suitable rescaling, the profile of the DoS tail depends only on dimensionality and 
separation from the support. In particular, in contrast to Lifshitz states in the Hermitian 
system, properties of the tail states are universal, independent of the nature of the random 
distribution. 

To illustrate the generality of our approach we will consider two examples which, 
according to the classification scheme defined below, belong to different fundamental 
symmetry. Firstly, following our discussion above, we will investigate the properties of 
tail states in a quantum system involving the propagation of a particle in a random scalar 
potential involving both real and imaginary components (|^). Secondly, we will investigate 
the nature of tails states in a quantum system where the particle is subject to a random 
real scalar potential and a random imaginary vector potential. 

As well as the Hatano- Nelson system where h = const, (|^) has been studied in the 
context of random classical dynamics involving classical diffusion in the background of a 
quenched random velocity field. Identifying I /2m with the classical diffusion constant Z?, 
and h/2m with a random velocity field v, a redefinition of the random potential W{r) 
obtains the classical operator 

H = -DV^ + (V • v) + V V + W{r) . 

In the quasi-classical approximation, (i.e. setting z as the complex energy, where 
r3ff [z] h) our goal will be to show that the long-range, low-energy properties of the 
random systems are described by a supersymmetric non-linear cr-model. At the level of 
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mean-field, the spectrum of both non-Hermitian operators will be shown to be bound to 
a region of support inside the complex plane. By taking into account non-perturbative 
symmetry breaking instanton field configurations of the action, we will reveal the existence 
of tail states which extend beyond the region of support predicted by the mean-field. We 
will determine the localization properties and profile of the DoS revealing an intrinsic 
universality of the present scheme. 

The paper is organized as follows: after introducing a general scheme to classify lin- 
ear non-Hermitian operators on the basis of their fundamental symmetries, in section ^ 
we formulate a statistical field theory to describe the spectral properties of the random 
imaginary scalar potential Hamiltonian. Within this approach, we show that the large 
scale properties of the spectral support are determined by the saddle-point or mean-field 
properties of the theory. Motivated by a parallel description of the disordered supercon- 
ducting system, in section 2.3 we show that optimal fluctuations of the random impurity 
potential induce localized tail states which extend beyond the region of support. Within 
this approach, we determine analytical expressions for the spatial extent of the localized 
states in the tail region, and determine the scaling of the complex DoS as a function of 
energy. Generalizing the scheme to the consideration of the imaginary vector potential 
Hamiltonian in section ^, we emphasize the universality of the present scheme. Results 
from both sections are compared with numerical simulations. Finally, in section ^ we 
conclude our discussion. 



2 Field Theory of non-Hermitian Operators 

The aim of this section is to prepare a statistical field theory of the weakly disordered non- 
Hermitian system. To be concrete, in the first instance, we will focus on the problem of the 
imaginary scalar potential (^. Later, we will see that the field theoretic scheme is easily 
modified to accommodate the constant vector potential perturbation. However, before 
embarking on this program, it is useful to first contrast general properties of complex 
linear non-Hermitian operators with those of their Hermitian counterparts. As well as 
helping to clarify the idiosyncrasies of the non-Hermitian system, we will expose a general 
scheme in which non-Hermitian operators can be classified according to their fundamental 
symmetries. 

2.1 Background: Symmetry Classification 

A general non-Hermitian operator H is specified by left and right cigcnfunctions, 

H\Ri) = {L,\H = {U\z^ , 

where {zi} denotes the set of complex eigenvalue. Although not a generic feature of non- 
Hermitian operators, throughout we will rely on the assumption that the eigenfunctions 
form a complete basis set. This is the case if there are no repeated eigenvalues, which will 
be true for a generic random operator such as those studied here. Furthermore we will 
take the eigenstates to be orthonormal, 

Y.\Rr){L,\=l. (3) 

i 

Now as with Hermitian operators, spectral (and localization) properties of the non- 



{Li\Rj) = Sij 
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Hcrmitian Hamiltonian can be expressed through the complex Green function 



9{z) = J- , 

z — H 

where z = x + iy denotes the complex argument. Inserting the resolution of identity (^), 
one obtains the spectral decomposition, 

giz) = Y,\R,)^{U\. (4) 

^ Z- Zi 

Using the analytical properties of the resolvent, it is straightforward to show that the 
density of complex eigenvalues 



.(z) = -^tr<5(z-i?) = -l^lim5: 



,2 



(|z-Z,P+,y2)2 



is expressible in terms of the complex Green function through the identity (Q): 

i^iz) = ^^2* ^rg{z) . 

The complex Green function g{z) is, therefore, a non-analytic function everywhere in the 
complex z-plane in which the corresponding eigenvalues density is non- vanishing. 

Previous studies have shown that techniques based on standard diagrammatic per- 
turbation theory account only for contributions to g{z) which are analytic in z [p] |. To 
account for all contributions, it is convenient to follow the now standard route |27, ^, 



of expressing the Green function through an auxiliary Hamiltonian which is explicitly 
Hermitian. This is achieved by constructing a matrix Hamiltonian, 7i, with the following 
2x2 block structure 



n 



H - z 

-z* 



In this representation, the Green function of the non-Hcrmitian operator is expressed as 
the off-diagonal element of the matrix Green function , 

g{z) = lini G21 (?7, z) , 

where G{r], z) = {irj — 7i)~^. 

This 'method of Hermitization' affords a convenient way of classifying the symmetries 
of a non-Hermitian Hamiltonian according to the fundamental symmetries of its Hermi- 
tian counterpart. The utility of this classification will become manifest presently in the 
construction of the statistical field theory of the non-Hermitian system. As an example, 
let us consider the Hamiltonian involving an imaginary scalar potential (|l]). Applied to 
this problem, the Hermitization procedure leads to the matrix Hamiltonian 

n = (^Cp - ^ + C7i + {y ~ V) a2 , 

where the Pauli matrices cr^ act in the 2x2 sub-space, and Cp = fp' /2m. 
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In fact, the matrix structure of 7i bares much in common with the Gor'kov or Bogohu- 
bov-de Gennes quasi-particle Hamiltonian of a weakly disordered superconductor [^sj . To 
foster this analogy let us implement the canonical transformation 

n^nr = T-^nr = (Cp-x + w)a3 + {y- V) (T2 , (5) 



where F = exp[— i7r(T2/4]. Associating the matrix structure with a 'particle/hole' space, 
the operator Tir can be identified as a Gor'kov Hamiltonian for a disordered supercon- 
ductor. Here x plays the role of the chemical potential, while Cp + W denotes the bare 
Hamiltonian, and {y — V) represents the (in this case real, random) order parameter. 
However, taking the analogy further, the energy argument in the Gor'kov Green function 
corresponds, in this case, to the infinitesimal irj of the matrix Green function Gj^^{rj,z) 
— evidently, we are interested in the zero energy quasi-particle states of the disordered, 
time- reversal invariant, superconductor with a random real order parameter. 

In its Hermitian form, the Hamiltonian above can be classified according to its fun- 



damental symmetries. According to the classification introduced by Zirnbauer |36 , the 
matrix construction places a general non-Hermitian Hamiltonian in the chiral unitary 
symmetry class denoted AIIl|^, a general element of which has the form 

Z 

Zt 

where Z is an arbitrary complex matrix. However, within this group, there are four sub- 
classes of higher symmetry: class CI (the time-reversal invariant superconductor) where 
matrices Z are complex symmetric; class BDI (chiral orthogonal) where the elements of 
Z are arbitrary and real; class CII where Z has an underlying symplectic structure; and 
class Dili where Z is complex antisymmetric. As will be clear from the forthcoming anal- 
ysis, the statistical field theory describing the non-Hermitian spectral properties has soft 
modes associated with each universality class. From our discussion above, the imaginary 
scalar potential is accommodated in the symmetry class CI. 

As the second example, let us consider the non-Hermitian Hamiltonian (^) involving 
the random imaginary vector potential. In this case, the Hermitization procedure leads 
to the matrix Hamiltonian 

«-(^— )'-(-^^)- 

As a result, we find that the system belongs to a different class depending on the value of 
the imaginary component of the energy argument y. Along the real axis (i.e. for y = 0) 
H is real and therefore belongs to the symmetry class BDI. Away from the real axis the 
Hamiltonian becomes complex and the symmetry is lowered to class AIII. Later we will 
discuss the physical manifestations of the discrete symmetries in the two limits. 

2.2 Generating Functional 

With this background, let us now turn to the construction of a field theory to describe 
statistical correlations of the non-Hermitian model. Motivated by the correspondence 
outlined above, the analysis will parallel previous investigations of the weakly disordered 



^ The notation adopted by Zirnbauer is motivated by that introduced by Cartan to classify the 10 -f 1 
symmetric spaces. 
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superconductor 



p8\ |39| | . The starting point is the generating functional for the single- 
particle Green function. Adopting throughout the convention H = c = 1, the latter is 
expressed as a field integral involving two independent four component supervector fields 
ip{r) and V'(r) with elements having both commuting/anti-commuting and 'particlc/hole' 
indices: 



D{ip, ip) exp 



dr 



Ip {i-q — l-djip+ / dr {%pj -f jif) 



Elements of the Green function can be generated from the source terms j'(r) and j(r). For 
example, from the generating function one can obtain the DoS according to the identity: 

dr 6^ 



v{z) 



—dz- lim / 

TT J 



Ld 6j^ir)Sjf{r) 



m 



(7) 



To be concrete, let us first consider the properties of the Hamiltonian involving the 
imaginary scalar potential. To exploit the analogy with the superconducting system, 
it is again convenient to implement the gauge transformation ijj i— > Tip whereupon the 
generating functional assumes the form 



Z[0]= / D{^,i^)exp 



i / dr 'tp{r) [irj — Hr) fpir) 



(8) 



where Tir is defined above (^). Now, as with the superconductor, the matrix Hamilto- 
nian ^ satisfies the particle/hole or 'charge conjugation' (cc) symmetry 

— (T2^rO'2 = • 

The latter is responsible for quantum interference effects which modify the long-range or 
low-energy properties of the average Green function. To accommodate these effects it is 
convenient to double the field space 

2ij [iT] - Tir) = {iV - ^r) ^ + V""^ [ir] - Upj V 

^ tp (^iri- Tir^ ip + {irj + (T27irCT2^ 'ip^ 

= 2* (iriaf - Hr) * , 

where, defining the Pauli matrix 173'^' which operates in the charge-conjugation (cc) space. 



1 



1 



As a result, the generating functional takes the form: 



Z[0] = / D{^, exp <^ i / dr 5'(r 



C,t.-W\a:i^{v~V)a2 *(r) 



With this definition, the fields ^ and are not independent but obey the symmetry 
relations 

\I/ = cr27'I'^ ^ = -*^cr27^ , (9) 

with 7 = E'bb ® (7^° — ii^FF ® o'2*^, where £^bb = diag(l, 0)bf and -Eff = diag(0, 1)bf are the 
projectors on the boson-boson and fermion-fermion sector respectively. This completes 
the construction of the generating functional as a functional field integral. To make further 
progress, we can draw on the intuition afforded by existing studies of the superconducting 
system. 
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2.2.1 Disorder Averaging 



As a first step towards the construction of an effective low-energy theory it is necessary 
to subject the generating functional to an ensemble average over the random impurity 
distribution. For this purpose we will take both the real and imaginary components of 
the random scalar potential to be Gaussian (5-correlated with zero mean, and correlation 

^ -<5'^(ri - r2) {V{r,)V{r,))y^-^S''{r,-r2), (10) 



{Wiri)W{r2))w = 



where v = (AL'')"^ ~ m(2ma;)'^''^^-'/^ represents the unperturbed DoS of the Hermitian 
Hamiltonian at energy x, while r and t„ denote, respectively, the mean scattering time 
of the real and imaginary components of the potential. In the following, we will suppose 
that the imaginary component of the scattering potential is weak, i.e. t„ 3> t. 

The ensemble averaging over the real random scalar potential W^(r) induces an inter- 
action of the fields which can be decoupled by the introduction of an 8 x 8 component 
Hubbard-Stratonovich field Q, 



(exp 



dr ^Wa-i^ 



iw 



DQ exp 



dr 



8t 



strQ^ 



1 

27 



Transforming the fields according to the symmetry (||), one finds that the supermatrix 
fields Q{r) must be subjected to the constraint: 

Q-ai7QV^^i- (11) 

Similarly, the interaction of the fields induced by the ensemble average over the ran- 
dom imaginary scalar potential V{r) can be decoupled by an 8 x 8 component Hubbard- 
Stratonovich field P(r), satisfying the symmetry constraint, P = ^P^^~^: 



(exp 



dr * 1/(72* 



dr 



^strp2 - -^^'Pcr2^' 
8t„ 2t„ 



Finally, an integration over the fields ^' obtains 



)w,v 



= DQ DPcxp 



— / dr str — 
8 ./ \ T 



P2 



1 



— - - / dr str(r|lng" 



where 



represents the supermatrix Green function. Further progress is possible only within the 
a saddle-point approximation. Following Ref. in the quasi-classical approximation 
X 3> 1/t ^ [1/T„,y], it is convenient to implement a two- level procedure in which we 
first formulate an intermediate energy scale action in which the terms in y and 1/t„ are 
treated as small symmetry breaking sources. Later, the influence of V (through P) on 
the complex DoS can be explored within a stationary analysis of this reduced action. 

2.2.2 Intermediate Energy Scale Action 

Neglecting the terms in y and P, a variation of the action with respect to Q obtains the 
saddle-point equation 



Q(r) = — eo(r,r) 



Qq ^ = ^770-3^ (g) 0-3 
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In the quasi-classical limit, xt ^ 1, taking Q to be homogeneous in space, this equation 
can be solved in the pole approximation with Q = ±1. Taking into account the analyt- 
ical properties of the average Green function, we identify the solution Qsp = (J^® Cg^. 
However, in the limit 77 0, the saddle-point solution is not unique, but expands to 
span a degenerate manifold of solutions generated by the homogeneous pseudo-unitary 
transformations Qs-p = Ta^ ® a^°T~^ compatible with the symmetry of Q, (p]). 

Fluctuations around the saddle-point can be classified into longitudinal and transverse 
modes according to whether or not they violate the non-linear constraint Q^(r) = 1. In 
the quasi-classical limit, the longitudinal fluctuations are rendered massive and do not 
contribute to the low-energy properties of the generating functional. On the other hand, 
transverse fluctuations of the fields become soft and must be accommodated. Taking into 
account soft fluctuations Q(r) = r(r)cr3 (g) a^'^T~^(r), and restoring the perturbations y 
and P, a gradient expansion of the action obtains the average generating functional 



= / DQ DPexp^— dr strP 



exp < — / dr str 



^2 



8t„ _ 

D{VQf + Ai (ir]az ® af' - iyai - ^cti ) Q 



where D = v^r/d is the diffusion constant. Here we have neglected terms higher order in 
yr and r/Tn- In the same approximation yr ^ 1, the supermatrix Green function takes 



the form 



^o(r,r') = -*Wd(|r-r'|)Q 



r + r 



2 

where fd{r) = T{d/2)[2/kFr]'^^^'^Jd/2~i{kFr)ex.p[-r/2£] is the Friedel function, and 
kp ^ \/2mx. 

Finally, integrating over the supermatrices P, we obtain the average generating func- 
tional 



where 



/ DQexp{~S[Q]} 



S[Q] = y dr str 



D{VQf - 4 (77^3 ® af - yai) Q + — [a^Qf 



Tri 



(12) 



This completes the construction of the intermediate energy scale action for the non- 
Hermitian system. Returning to the analogy with the superconductor, we note that the 
form of the action (|l^) coincides with that of a disordered superconductor subject to a 
random real order parameter with an average magnitude y and variance set by 1/t„. To 
investigate the influence of the imaginary potential on the complex DoS, it is necessary 
to subject this action to a further stationary analysis and obtain the low-energy action. 

2.2.3 Low-Energy Scale Action 



Although the soft mode action (12) is stabilized by the quasi-classical parameter xt, the 



terms in y and 1/t„ render the majority of field fluctuations massive. To obtain the 
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low-energy saddle-point, let us again vary the action with respect to Q subject to the 
non- linear constraint, = I: 



DViQVQ) ~ [rycTa ® of - ya^,Q\ -f -^[ctiQcti, Q] = . 



Applying the Ansatz Q — cosOas (E) erf -|-sin6'cri, where 9 = diag(i6'BB, ^ff)b,f, the saddle- 
point equation assumes the form 



7] sm t) + y cos ( 



— sin2^==0, (13) 

Tn 



a result reminiscent of the Abrikosov-Gor'kov mean-field equation of the disordered su- 
perconductor subject to a pair-breaking perturbation pO| . Taking the solution to be 
symmetric and homogeneous in space, one obtains 

- arcsin(yT„) |y|T„ < 1 

(14) 

--sign(?/) |y|r„ > 1 . 

Once again, as rj — > 0, this solution is not unique but expands to span an en- 
tire manifold parameterized by transformations Q = TQ^^T^^ , where [r, cti] = and 
T — 7(T^^)^7^ — these represent the soft fluctuations corresponding to the time-reversal 
invariant superconductor symmetry class CI. Substituted back into the action, these fluc- 
tuations are described by the corresponding low-energy action 

SosiQ] j dvD{l- y^rl) str(VQ)2 , 

where Q = Tct3 ® af 'T~^ . This completes our construction of the field theory describing 
spectral correlations in the disordered non-Hermitian system. In the following section we 
will use these results to explore the support of the complex DoS. 

2.3 Complex Density of States 

Making use of Eq. (Q) , the complex DoS is obtained from the generating function as 

{v{z,v))^y = ~'-vd,^, lim(str [(I - ai) ® E}}, ® afQ{v)\)Q , (15) 

where {■ ■ ■)q = /q2^j DQ . . . e"'^!'^! where S[Q] represents the intermediate energy scale 
action (|l2|). In the homogeneous saddle-point (i.e. mean- field) approximation (pH), the 
corresponding DoS takes the form 



\y\rn < 1 

{i'{z))w,v — ivdz* lirn sin^MF — (16) 

\y\Tn > 1 . 



As expected, at the level of the mean-field, the field theory predicts a migration of the 
DoS off the real line and into the complex plane. The density of complex eigenvalues 
is constant over the region of support set by the effective scattering rate of the non- 
Hermitian potential. Reassuringly, the expression for the DoS satisfies the sum rule 
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Figure 2: Complex DoS in the quasi-classical limit xt ^ 1 for the two-dimensional scalar 
potential problem, (a) In the mean-field approximation, the DoS is equal to VTn/2 inside the 
region |j/|t„ < 1 of the complex plane z = x + iy and zero outside; (b) {i/{z)}w,v versus y for a 
fixed X and for two values of the scattering time, t„ > r^: while in the mean- field approximation 
th e D oS shows sharp edges, rare realizations of the random potential give rise to tail states 
(§11). 



J dy{i^{z))w,v = V. In the two-dimensional case, the result is illustrated qualitatively in 
Fig. g. As the scattering time t„ does not depend on the real part of the energy x, the 
width of the mean- field spectrum boundary is constant. 

To complete the analysis, it is necessary to take into account the fluctuations of the 
fields around the saddle-point. As discussed above, these fluctuations divide into a set 
which are massive (on the scale of y) and a massless set. Interestingly, one finds that the 
latter commute with the source and do not contribute to the DoS! 

Indeed, the general profile of the support is born out by numerical simulation. Figure |l| 
shows the amalgamation of data for the two dimensional random system generated from 
100 realizations of the complex random scalar potential Hamiltonian. However, instead 
of a sharp cut-off in the support, the data clearly indicate a soft edge with eigenvalues 
occupying the region outside the mean- field prediction (|l^). How can this observation be 
reconciled with the results of the field theory? Surprisingly, taking into account massive 
fluctuations of the action perturbatively, the mean- field prediction remains intact. In fact, 
the tail states that decorate the edge of the support are generated by optimal fluctuations 
of the real random impurity potential and are reflected in non-perturbative instanton 
configurations of the action. 



2.4 Tail States and Instantons 

To identify corrections to the mean-field DoS ( p^ ) we can draw on the intuition afforded 
by recent studies of sub-gap state formation due to optimal fluctuations in the supercon- 
ducting system While the homogeneous solution of the mean-field equation gives 
rise to a hard edge of the DoS support, we will see that inhomogeneous symmetry broken 
field configurations reflect the influence of rare realizations or 'optimal fluctuations' of 
the random scalar potentials which soften the edge (Fig. |2|). Here, for simplicity, let us 
focus on the quasi one-dimensional case, generalizing the discussion to the d-dimensional 



system in section 2.4.2 
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Figure 3: Integration contours for boson-boson and fermion-fermion fields in the complex plane 
9' = —i9 sign{y). The two degenerate bounce solutions are qualitatively shown, together with 
the mean- field starting point (/)mf and the inversion point (j)inv- 



2.4.1 Quasi One-Dimensional Geometry 

To explore the influence of inhomogcneous field configurations of the action, let us revisit 
the mean-field equation of motion ([T^). Operationally, it is convenient to deal not with 
the saddle-point equation ( |l3|) itself, but rather its first integral, D (V 6)^ ~ yV (6) — const, 
where 



V{§) 



■ cos 29 



yTn 



denotes the effective complex 'potential'. 

To identify the inhomogeneous instanton solution let us first note that, outside the 
region of support, |y|r„ > 1, the homogeneous mean- field solution takes the form ^mf = 
—tt/2 sign(?/). Taking into account the condition that the solution should coincide with 9uf 
at infinity, one can identify a 'bounce' solution parameterized hy — [— 7r/2 -I- ic/)] sign(2/), 
with (f) real, involving the real potential ^4(0) = V{[—t:/2 + i0]sign(y)) with endpoint 

0inv such that VR(0inv) = K(0mf)- 

Now integration over the angles is constrained to certain contours . Is the bounce 
solution accessible to both? As usual, the contour of integration over the boson-boson 
field 9bb includes the entire real axis, while for the fermion-fermion field i0pp runs along 
the imaginary axis from to in. With a smooth deformation of the integration contours, 
the mean- field saddle-point is accessible to both the angles [|7[ |3^. By contrast, the 
particular bounce solution and the mean-field solution can be reached simultaneously by 
a smooth deformation of the integration contour only for the boson-boson field 0bb (see 
Fig. The particular bounce solution breaks supersymmetry. 

Applying the parameterization for the 'bounce' 

iOBB{r) =12^ sign(y) 6'ff = 9mf , 

the first integral takes the form 



(V,/j0)2 + K(0) = K(0mf) 



= -4cosh(/)-H/t"icosh20 , (17) 
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Figure 4: Rescaled potential Vr{(I)) — Vr{0) versus the variable (j) for k — 2.5, k — 1.6 and k = 1. 
where ^ = {D/\y\y^^ and 

K = r„|y|. (18) 

Here (^mf — for |j/|t„ > 1, while (f>MF is purely imaginary for |i/|r„ < 1. Typical shapes 
of the potential Vn{4>) for different values of k are shown in Fig. |[ A bounce solution 
exists only for values of n bigger than one, while for n < 1 the unique solution is the 
homog eneous one, cp — (Pmf- 

For the quasi one-dimensional system it is possible to derive analytic expressions for 
both the real instanton action 5inst — '2T^i^\y\^Stf,{K) 



S4,{k) = / d0v/K(O)-K(<^) , (19) 







where i^inv = | arccosh(2K — 1)| is the endpoint of the motion, and the bounce solution 
^(r). In particular, for S^{k) we have (Fig. ^ 



^^(k) = %/2^|-2 ^^^ - + 2arctanV^7^| . (20) 

At the same time, imposing the boundary conditions (^(r — > ±oo) = 0, from the first 
integral (^7|) one obtains the bounce solution 

, , K2-)-2e~2|'-|/'-oW^Ar^(3K-2)-|^e~4|r|/ro(«)(-^_2) 

cosli 0(7*1 — ■ 1 

K2^_2e-2k|/r„(K)^/^7ZT-(2-K)+e-4kl/'-o(«)(,i-l) ' 

where 

r dt 1 
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Figure 5: Exact (solid line) and approximate (dashed line) bounce solution (j)(r) versus 
r / y' Dt„/2 in the quasi one-dimensional geometry for two different values of the inverse of the 
imaginary 'impurity' concentration, k = 1.6 and re = 1.05. The interval [— ro(/v), ro(ft)] is explic- 
itly plotted for the exact solutions. 



defines the extent of the 'droplet' (see Fig.|5|). In particular, we note that, one approaching 
the band edge (k — ^ 1+), ro(K) oo. 

Now close to the band edge {k 1+), the maximum of the potential Vr((/)) converges 
on the two minima allowing for the development of a perturbative expansion in ~ 0mf = 
0. Setting /? = (4 - k)/6k ~ 1/2, 

K(0) ^ K(0) - 2^02 ^ ^^4 ^ 0(^6) (22) 
In this approximation the action and the bounce solution are given respectively by: 



S,{.) ^(«-l)3/^ Hr) - 4V2/" 



K>i 3 K>i V K e-2k|/ro(K) _|_ 4^ ■ 

This completes the estimate of the contribution to the action from a single bounce 
configuration. However, to complete the analysis it is necessary to explore the infiuence 
of fluctuations around the instanton solution. Here we sketch the important aspects of 
the analysis referring for details to the parallel discussion by Ref. |4l| in the context of 
the disordered superconductor. Generally, field fiuctuations around the instanton solution 
can be separated into 'radial' and 'angular' contributions. The former involve fluctuations 
of the diagonal elements 9, while the latter describe rotations including those Grassmann 
transformations which mix the bf sector. Both classes of fluctuations play an important 
role. 

Dealing first with the angular fiuctuations, supersymmetry breaking of the bounce 
is accompanied by the appearance of a Grassmann zero mode separated by an energy 
gap from higher excitations. This Goldstone mode restores the global supersymmetry of 
the theory. Crucially, this mode ensures that the saddle-point respects the normalization 
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Figure 6: Exact action 5^ (solid line) versus k, in the quasi one-dimensional geometry together 
with the first term of the expansion in (k — 1) (dashed line). The action vanishes as k ^ . In 
the inset the rescaled DoS {v(£))w,v jv ~ exp[— Sinst] versus j/Tn, for -kv^/dJt^ = 2, is plotted. 
For values of Ksign(y) near the critical values ±1, the saddle-point analysis becomes unreliable. 
Here, fluctuations lead to a smooth interpolation between the DoS in the tail region, and that 
coming from the bulk mean-field. 



condition (2^[0])m/,\/ — 1. Associated vi^ith radial fluctuations around the bounce, there 
exists a zero mode due to translational invariance of the solution, and a negative energy 
mode (c.f. Ref. H). 

Combining these contributions, one obtains the following expression for the local com- 
plex DoS in the tail region, 

{iy{z,r))w,v ii'dz- [cosh0(r) - cosh^Mp] |xo(?-)Pe-'^'°=* , (23) 

K>1 

where Xo('') represents the eigenfunction for the Grassmann zero mode and Siast denotes 
the instanton action (pO|). Thus, to exponential accuracy, the complex local DoS in the 
tail region becomes non-zero only in the vicinity of the bounce. Mechanisms of quantum 
interference due to optimal potential fluctuations in the non-Hermitian system conspire to 
localize states in the tail regions on a length scale ro{K) greatly in excess of the wavelength 
A = l/kp. 



2.4.2 Generalization to Dimensions 1 < < 4 



Close to the band edge (k > 1), a generalization of the quasi one-dimensional results to 
higher dimensions can be developed by dimensional analysis. From the expansion 
for the potential Vr ((/>), the bounce conflguration is shown to have the scaling form 
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Figure 7: Plot of the (logarithm of the) complex DoS versus Q[z] = y for a fixed value of 5R[z] = x 
taken from the numerical simulation of Fig. |l| The exponential fit of the data is shown as solid 
curve. 

which, in turn, implies the instanton action (p^): 

S,^s, 4vra,(Ar„)('*-2)/2//2 [2(r„|y| - . (24) 

Here g = vDL'^^^ and A = {i'L'^)~^ represent respectively the dimensionless conductance 
and average level spacing of the Hermitian system (i.e. with = 0), and ad is a numerical 
constant (ao = 1/4 and ai = 1/3)0. Later, in discussing the universality of the results, 
we will return to consider the zero-dimensional situation. 

Taken together with Eq. (^3|), this result shows that the complex DoS in the tail 
region varies exponentially with the separation from the mean-field edge. The power law 
dependence of the exponent depends on dimensionality, presenting a linear dependence 
on (t„|j/| — 1) in the two-dimensional system. In this case, moreover, the only dependence 
in (|4|) on the real part of the energy a: occurs through the dimensionless conductance, 
g = vD ~ I a; I, so that the width of the exponential tail scales as \/\x\, as qualitatively 
appears in Fig. |^. A comparison with the numerical simulation is shown in Fig. |7[ The 
results show a good quantitative agreement with the theory for the tail of the DoS. Note 
that close to (but within) the edge of the support of the spectrum, the saddle-point 
approximation used in estimating the mean-field DoS becomes unreliable. As shown by 
the zero-dimensional analysis below, an honest treatment of fluctuations in this region 
provides a smooth interpolation between the bulk and tail states. 

3 Random Imaginary Vector Potential 

To conclude our discussion of tail states in the non-Hermitian system, we turn now to 
an analysis of the random imaginary vector potential Hamiltonian (H). Once again, our 



2 Specifically = /du[(Vu/)2 + /^(u) - /4(u)]. 
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approach will be to formulate an effective field theory from which the spectral proper- 



ties can be determined. Although, as shown in section 2.1, the random vector potential 



Hamiltonian belongs to a different fundamental symmetry from the random scalar poten- 
tial Hamiltonian, we will show below that, with some qualification, the DoS in the vicinity 
of the support exhibits the same universal scaling form. 

To be concrete, let us consider a d-dimensional Hamiltonian (|^) involving an imagi- 
nary vector potential field h(r) together with a real random scalar potential W{r). Once 
again, let us assume that the scalar potential W{r) is drawn from a Gaussian (5-correlated 
impurity distribution (|lo|). Similarly, we will be primarily concerned with a system ex- 
hibiting a random vector potential drawn from a Gaussian distribution with zero mean 
(h(r));i = and correlation 



(/ii(q)/i,(q')>/.-(2^)''<5(q + q')/(|q|) 



(25) 



Here, for reasons that will become clear below, we have introduced an 'envelope' /(|q|) 
which limits correlations of the fields at both microscopic and macroscopic length scales. 
Setting 



h(r) = V^{v) + V A x(r) 



(26) 



7i is identified as the strength of fluctuations of the incompressible component of the 
vector field, x(r), while 72 controls the irrotational part, i^(r). In the following, it will 
also be useful to contrast the behaviour of the random vector potential system with that 
of a constant vector field h(r) = const = ho — the Hatano-Nelson system H. 



3.1 Background 

Before turning to the formal analysis of the statistical field theory, which begin with some 
general considerations which identify certain idiosyncrasies of the constant vector poten- 
tial system and how they impact upon the existence of tail states in the random vector 
potential system. Previous investigations of the former have revealed unusual localiza- 
tion properties of the random Hamiltonian which contrast with those of its Hermitian 
(i.e. ho — 0) counterpart. Specifically, it was shown in Ref. |^ that, for the strictly one- 
dimensional system subject to periodic boundary conditions at infinity, when the field, |ho|, 
exceeds a critical value he corresponding to the (energy dependent) localization length 
Cioc = 1/^c of the Hermitian system, there is a transition to a dclocalizcd phase. Since 
then, this result has been generalized to the quasi one-dimensional system by Kolesnikov 
and Efetov ||l^ using techniques which parallel those discussed here. 

Now, following Ref. 1^, a simple argument can be established to describe qualitatively 
the origin of the transition. For a constant imaginary vector potential, ho, a similar- 
ity transformation of the left and right eigenfunctions 4>f{r) t-^ (r) exp(— ho • r) and 
4>f{r) 1-^ exp(ho • r)(j)f-{r) removes from the Hamiltonian its dependence on ho. Being 
now Hermitian, the eigenvalues of the Hamiltonian must therefore be real and specified 
by those of the unperturbed random system. However, if periodic boundary conditions 
are imposed, the similarity transformation must be applied with caution: if the eigen- 
functions of the Hermitian system are localized on a length scale ^loc < |ho| the similarity 
transformation is incompatible with the periodic boundary conditions and can not be 
imposed. Here there is a transition from a localized phase where the eigenvalues are real 
to a delocalized phase where the eigenvalues migrate into the complex plane. 
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Figure 8: A typical field configuration of the random vector field. Over the interval £h, the flux 
lines of h are approximatively constant, preventing the existence of tail state. 

Now what does this phenomenology tell us about the existence of tail states? Using the 
same arguments, it is easy to deduce that optimal fluctuations of the impurity potential 
can not induce localized tail states in the constant vector potential system: suppose that 
a tail state appeared at some complex energy outside the band of bulk states. If the state 
is localized, it must be insensitive to the boundary conditions at infinity and, therefore, its 
dependence on the constant field can be removed by a similarity transformation. As such, 
it must therefore be contained within the Hermitian theory, and it's eigenvalue must be 
real. This contradiction strictly eliminates tail states from the infinite constant imaginary 
vector potential system. (Of course, the Hermitian system can and will exhibit low-energy 
Lifshitz band tail states.) 

With this background, let us now turn to consider the random vector potential sys- 
tem. Once again, for a random field h(r) drawn from the Gaussian distribution (p5|), it is 
straightforward to confirm that the similarity transformation (f>f{r) i-^ (/^^(r) exp[— (^(r)] 
and (pfij) ^ exp[(p(r)]0f (r) removes from the Hamiltonian the dependence on the incom- 
pressible component, ip{v). Therefore, when subject to any purely irrotational vector field 
h(r) = V(^(r) (with zero average), the eigenvalues remain real and coincide with those of 
the Hermitian system. Conversely, the incompressible field distribution h(r) = V A x(r) 
can not be removed by similarity transformation. Such field configurations necessarily 
generate states with complex eigenvalues. 

Now, the random distribution ( p5[ ) with /(|q|) = 1 imposes long-ranged, i.e. power- 
law, correlations in the incompressible field configuration. As well as admitting superdif- 
fusion processes ^ which dominate the relaxation to equilibrium in the classical 

system, these long-range correlations impact upon the probability of finding tail states in 
the non- Hermitian system: any random configuration of the incompressible field will in- 
volve field lines correlated over arbitrarily long distances. These field lines advect particles 
in the classical system and are responsible for the characteristic superdiffusion properties 
of the system. Now in such a background, one can identify regions where the field lines are 
oriented over long distances, say Ih (see Fig. ||). Within this region, the nucleation of tail 
states with a localization length ro < is prohibited by the same mechanism described 
above for the Hatano-Nelson system. Clearly, the phenomenology of tail state formation 
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in the superdifFusive system is subtle being sensitive to the arrangement of both the scalar 
and vector field. 

Therefore, since we are interested here in exposing the general principles behind tail 
state formation in the non-Hermitian system, to keep our discussion simple, we will limit 
our considerations to distributions where the correlations of the field lines are limited to 
a length scale £h much smaller than the localization length of the incipient tail states, tq, 
i.e. introducing a microscopic cut-off a, we set 



/(|q|) 



1 1/4 < |q| < 1/a 
otherwise , 



where £h represents the maximum length scale of field correlations. With these considera- 
tions in mind, let us now turn to the development of the field theory of the non-Hermitian 
system. 



3.2 Field Theory 

Following the method outlined for the imaginary scalar potential system, our starting 
point is the gauge transformed 'Hermitized' matrix Hamiltonian (^. Now, under the 
'charge-conjugation' operation, the Hamiltonian transforms as 



(72 [Hr(h)J (72 = -Wr(-h) . 

The statistical properties of the non-Hermitian system can again be obtained from the 
generating function (||) by forming an average over the real random potential VF(r) and 
introducing a slow field decoupling. In doing so, one obtains the generating functional 



{Z[0])w = J DQexp 



- /drstrg^-- 



dr str(r|lng"Mr) 



(27) 



where 



"1 — i-qa-i ® (73'' 



1 

2m 



[p - ih(r)cri 



2r 



Q(r) 



denotes the supermatrix Green function. As with the scalar potential, further progress 
is possible only within a saddle-point approximation. However, following our discussion 
above, it is now necessary to exercise some caution. 

In the Hermitian system with a real vector potential A = ih, the conventional ap- 



proach 1 35 involves subjecting the action to a saddle-point approximation in the absence 
of the field. The saddle-point manifold = I is unperturbed by the real vector potential 
and its effect on the low-energy properties of the system can be accommodated at the 
level of the gradient expansion. However, as emphasized in Ref. in the presence of 
a constant imaginary vector potential, h(r) = ho, depending on the magnitude of ho, it 
is possible to identify separate saddle-points of the total effective action (|27|). Inside the 
localized phase (i.e. |ho| < /ic), the dependence on hp can (and must) be removed by 
a rotation of Q — the counterpart of the similarity transformation. The resulting the- 
ory reflects that of the Hermitian model and the eigenvalues condense onto the real line. 
Conversely, in the delocalized phase (i.e. |ho| > he), the similarity transformation on Q 
is inadmissible. Here the dependence of the action on ho must be developed explicitly, 
reflecting the fact that the eigenvalues acquire an imaginary component. Similarly, for a 
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Figure 9: Complex DoS for the imaginary vector potential problem in the quasi-classical limit 
XT ^ 1 for the two-dimensional system, (a) In the mean-field approximation the DoS is equal 
to i'/{4:D'yi) inside the region \y\ < 2D71 = 2T'yi\x\/m of the complex plane z = x + iy and 
zero otherwise; (b) {u{z))w,h versus y for two values of 71 < 71. The tails of the DoS for 
\y\ > 2T'yi\x\/m are shown schematically. 



random vector potential, it is necessary to remove the irrotational component of the field 
h (i.e. V(^(r)) explicitly at the level of the saddle-point by subjecting Q to a similarity 
transformation 

Q(r) ^ e-'^('"''"i®'""Q(r)e'^('")''^®''" . 

— As expected, the same transformation leaves the source term for the DoS (|l5|) unper- 
turbed. 

Without the irrotational component of the field, the saddle-point of the total action 
can be analyzed in the usual way, following a two-stage procedure. As before, taking into 
account the quasi-classical parameter xr 3> 1, one finds that the low-energy fluctuations 
are confined to the manifold Q^{y) = I. Then, subjecting the action to a gradient expan- 
sion keeping yr and ^ £^^, one obtains the non- linear ci- model action ^, |8, ^ 

S[Q] = J dr stv [D{VQf - 4 (77^x3 ® of - yo^) q] , (28) 

where V = V — h[ , cri (g) af'\ represents the covariant derivative, and h = V A x reflects 
the incompressible coirrponent of the field. 

Now, although the structure of the action compares with that of the imaginary scalar 
potential field, it exhibits important differences which substantially influences the be- 
haviour of the system: specifically, in addition to the 'diamagnetic' contribution, i.e. the 
term proportional to h^, the action presents a paramagnetic contribution, linear in h, 
which couples to the gradient of Q. Although this term does not affect the homogeneous 
saddle-point solution, it is responsible for the destabilisation of the tail states in the sys- 
tem with long-range correlations. To see this explicitly, let us return to the consideration 
of the constant vector potential system. Previously, in the qualitative discussion above, it 
was argued that such a system does not accommodate tail states. How is this phenomenon 
reflected in the field theory? 

In fact, for h(r) = ho > h^, the action is still given by the non-linear cr-model (p^). 
However, in view of the paramagnetic term, to explore the saddle-point structure, it is 
necessary to adopt a more general parameterization, Q = (J\ cos B -f sin Be^"^ 0-3 ® af , 
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where p = diag(pBB, Pff)b,f- In doing so, two coupled equations for and p are obtained: 
i?V20 + :^sin2e[(Vp)2 + 4|hop+4ho- Vp] +2j/sin^ = (29) 



V • 



sin^ (V/5 + 2ho) 



= 



To understand the nature of these equations, let us focus on the quasi one-dimensional 
system. In this case the second equation can be integrated and substituted in the first 
one, giving respectively 

„ COS "* 

DV^e + Dj^ — + 2y sin 6* = Vp = -2ho 



sm^ sir? ( 



where j is a constant fixed by the boundary conditions. From these equations one can 
identify a homogeneous supersymmetric solution, V/5 = and 

a j~yTn \y\Tn < 1 

cos6Imf = < . , , 11^1 (30) 
l^-sign(y) |y|T„ > 1 , 

where l/r„ — 2Z3|hop. Again, making use of Eq. (^, the homogeneous solution translates 
to a complex DoS (^6|), which is flat and non- vanishing only in the interval |?/|t„ < 1. 
Now, in the imaginary scalar potential problem, an instanton configuration of the saddle- 
point equation was signalled by the existence of tail states. What happens in the present 
case? Since, at infinity, the instanton solution is constrained to be the homogeneous one, 
the constant j is determined by the homogeneous mean-field solution, j = 2ho sm C7^jp . 



Therefore, using (30), one can see that outside the bulk of the spectral support j vanishes. 
Therefore, the equation of motion for 9 does not admit the existence of an instanton. As 
expected from the qualitative discussion above, tails are therefore excluded from the 
constant imaginary vector potential problem. 

Similarly, in the random vector potential system with long-ranged correlations (i.e. 
for /(|q|) = 1), we can expect the paramagnetic term to present a similar role. Indeed, 
subjecting the generating functional to an average over the random distribution (^), 
one obtains a non-local interaction of the fields. In the present context, such terms have 
been shown by Taras-Semchuk and Efetov p6| to reproduce the known renormalisation 
properties of the superdiffusive system. However, in the present case, we have limited 
our considerations to correlations of the vector field which extend over a maximum range 
ih- Then, taking the relevant field configurations of Q (i.e. those involving the instanton 
solution) to be long-ranged on a scale tq greatly in excess of £h, the vector field entering 
the paramagnetic term can be replaced by its spatial average, i.e. the paramagnetic term 
can be eliminated. By contrast, the diamagnetic contribution survives spatial averaging. 
The total action assumes the form 



S[Q] = J dr str 



DiVQf - 4 (?7fT3 ® - yai) Q + — (Qdi ® of f 



(31) 



where, setting 71 = jiL'^ J dq/(27r)'^ /(|q|), 

— -2(d-l)Z?7i . 

Now, as with the imaginary scalar potential, the profile of the bulk DoS and tails 
states can be obtained by subjecting the action to a saddle-point analysis. Varying the 
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Figure 10: Complex eigenvalues for the for several realizations of the non-Hermitian vector 
potential Hamiltonian, where the potential W is drawn from a Gaussian (5-correlated random 
impurity distribution. Here we have included a lattice of 29 x 29 q-points in the momentum 
space. Notice that the constraints on the numerical simulation limit the relevant data to the 
interval x = < 200. 



action with respect to Q, subject to the non-linear constraint, = I, one obtains the 
homogeneous saddle-point equation, 

DW{QWQ) - [770-3 ® cr^ - y<^i,Q] + 7r-[^i ® ct^'^Qcti (g) 0-3'% Q] = , 

Then, applying the Ansatz Q = coaOa^ ® ctS*^ + sinOai, where 6 = diag(i0BBi ^ff)b,fi the 



saddle-point equation coincides with Eq. (y_3|) allowing the results of sections p.3| and 2.4.2 
to be imported. As a result, we can immediately deduce the homogeneous mean-field so- 
lution (Q) as well as the inhomogeneous instanton solution (|l3|) of the saddle-point equa- 
tion. Thus, from the homogeneous mean-field solution, one obtains the expression ( ]l6| ) 
for the DoS, i.e. the complex DoS is flat and non-vanishing over the interval 

Tn a m 

i.e., in the two-dimensional system, the support for the DoS occupies a wedge of the 
complex plane with a width that scales in proportion to 7i|a;| (see Fig. |^). Furthermore, 
making use of Eqs. (123) and (p4h, one obtains complex DoS in the tail region. 



{U{z))w,h ~ CXp {-4™d(AT„)('^-2)/2^^<i/2 J2(^„|y| _ l)](4-d)/2| _ 

where, as usual, g = vDL'^~^ and A = {yL'^)~^ represent respectively the dimensionless 
conductance and average level spacing of the Hermitian system. 

In fact, the behaviour of the vector potential system differs from the scalar potential 
system only in the symmetry of the soft fluctuations around the saddle-point solution. 
As usual, the homogeneous saddle-point solution Q^f is not unique but is spanned by an 
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Figure 11: (Logarithm of the) complex DoS versus Q[z] — y for a fixed value of = x taken 
from the numerical simulation of Fig. ^ (triangles). A second data set (circles) is shown for a 
value of 7i four times as large. Notice that, as expected from the theory, the change in 71 is 
reflected in the slope of the exponent. The exponential fits of the two data are shown as solid 
curves. 

entire manifold of solutions Q = TQmfT^^ where for y 7^ 0, [T, cti] = = [T, cri (g) erg''] 
with T = 7(T^i^)^7^ (Chiral symmetry class AIII in the classification of Ref. while, 
for y = 0, [T,ai (E) cr^'^] = (Chiral symmetry class BDI). Now while the soft (class 
AIII) fluctuations decouple from the DoS source (c.f. the scalar potential Hamiltonian), 
the soft (class BDI) fluctuations couple. As a result, the complex DoS becomes strongly 
suppressed near y = while, on the real axis, a finite density of eigenvalues accumulates 
— see below. 

The amalgamation of data for the eigenvalues of ca. 100 realizations of the imaginary 
vector potential Hamiltonian is shown in Fig. ^ The locus of the edge of the support 
predicted by the mean-field theory is in good agreement with the numerics, while the 
depletion of eigenvalues from the interval near y = and the singular density at y = 
are clearly resolved. In particular, a comparison of the complex DOS with the numerical 
simulation shown in Fig. |^ shows a good agreement with theory. 

3.3 Zero-dimensional Limit 

Finally, to conclude our discussion, let us turn to consider the properties of the imaginary 
vector potential Hamiltonian when the system enters the zero-dimensional limit. Our 
analysis above shows that, outside the region of support, the spectrum is characterized 
by tail states which are localized on a length scale of ro{K), (pl|). When the size of the 
droplet region becomes in excess of the system size (inevitable in a finite system as one 
approaches the edge of the support, k 1+), the system enters a zero-dimensional regime 
where the action is dominated by the zero spatial mode of Q. Here the action collapses 
onto the zero-dimensional theory analyzed in detail by Efetov Q . 

In this limit, we can proceed in two ways: firstly, following the discussion above. 
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we can apply the saddle-point analysis seeking a symmetry broken field configuration. 
However, in this case, the relative of the bounce solution is now a stationary point of the 
zero-dimensional potential V^ip). Alternatively, following Efetov Q one can undertake an 
exact evaluation of the zero-dimensional cr-model. A pursuit of the second route obtains 
the exact formula for the complex DoS. The latter can be separated into two contributions 
iy{z) — Vr{z) + Vc{z), where 

/ 2 2 TTi^ / "V \ / 2 2 

Vr{x,y) = v5{y) dt e""^ * iycix,y) = — $ [^j t sinh(xi)e-° ' , 



where 



2TT\y\ 

X = — ^ a 



2 



with — 2 dt e~* /^/^^ and A is the level spacing of the Hermitian system. The 
presence of the anomalous part Vrix, y) attests that a finite fraction of eigenvalues remains 
real. Otherwise the component Vc{x,y) describes the distribution of complex eigenvalues 
and its form is shown in Fig. |l2|. In the region x ^ 2a the DoS coincides with the 
constant mean- field value ( p^ for x < 2a^, while outside this interval one has the following 
asymptotic expansion: 

x>2a2 Ax V 4a"^ / 



With the exact result at hand, let us compare Eq. (|3^) with the result of the saddle- 
point analysis from the symmetry broken field configuration. In particular, in the zero- 
dimensional case, the saddle-point analysis leads to the 'instanton' action 

^inst = ^[K(0)-yK(^min)] , 
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where 0min is the minhxium of the potential Vr ((/)) (Fig. Evaluatmg this minimum, 
one obtains 

2 

From the expression for the DoS (p3|), we can conclude that, in the zero-dimensional case, 
the exponential dependence of the DoS coincides with the exact result obtained in Ref . |^ . 

4 Conclusions 

To conclude, we have implemented a field theoretic scheme to explore the structure of 
the DoS close to the edge of the support of two linear non-Hermitian operators describ- 
ing a quantum particle subject to a random imaginary scalar and an imaginary vector 
potential. In doing so, we have provided a general scheme for the symmetry classifica- 
tion of non-Hermitian operators. The field theoretic approach is easily generalized to the 
consideration of higher point spectral correlations of the fields. 

In the quasi-classical limit, where the real part of the energy is in excess of any other 
energy scale, the tails are dominated by 'optimal configurations' of the real random scalar 
potential. In contrast to band tail states in semi-conductors, these states are quasi- 
classical in nature being localized on the length scale ^ = {D/\y\y^^ '3> I- As such the 
profile of the DoS and their dimension is universal depending on just a few material 
parameters and independent of the nature of the impurity distribution. In the particular 
case of the constant imaginary vector potential, we have argued that tail states of the 
system are prohibited by the delocalization mechanism of Hatano and Nelson. 
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